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EFFECT OF RADIOMETRIC ERRORS ON ACCURACY OF TEMPERATURE- 
PROFILE MEASUREMENT BY SPECTRAL SCANNING USING 
ABSORPTION-EMISSION PYROMETRY 
by Donald R. Buchele 
Lewis Research Center 

SUMMARY 

The spectral -scanning method may be used to determine the temperature profile of 
a jet- or rocket-engine exhaust stream by measurements of gas radiation and trans- 
mittance at two or more wavelengths. A single, fixed line of sight is used, using immo- 
bile radiators and radiometers outside of the gas stream, and there is no interference 
with the flow. In the spectral -scanning method, the wavelengths of measurement are 
chosen to give greater weight to the radiation from particular portions of the optical 
path. Then, if a prior assumption is made that the temperature profile can be charac- 
terized by j independent parameters, these parameters can be determined from mea- 
surements at j wavelengths. A given radiometric error causes an error in computed 
temperatures . The ratio between temperature error and radiometric error depends on 
profile shape, path length, temperature level, the strength of line absorption, and the 
absorption coefficient and its temperature dependency. These influence the choice of 
wavelengths for any given gas. 

In the error analysis the fluxes are computed for an assumed profile shape. Differ- 
entiation then relates the radiometric error to the temperature error in terms of the 
more important experimental parameters of wavelength, absorption coefficient, and its 
temperature coefficient. Conditions for minimum temperature error are derived. Nu- 
merical results are presented for a two-wavelength measurement on a family of profiles 
that may be expected in a practical case of hydrogen-oxygen combustion. Under favor- 
able conditions, the fractional error in temperature approximates the fractional error 
in radiant -flux measurement. 



INTRODUCTION 


The spectral scanning technique in the absorption-emission method of gas pyrom- 
etry (ref. 1) provides a means of determining gas temperature along a path through the 
gas by using radiation sources and receivers located only at the ends of the path. In the 
case of a jet- or rocket-engine exhaust stream, the instrumentation may therefore be 
outside the stream. 

“ Any one measurement of radiant flux at the receiver represents the integral of 
radiant flux from every element of the path length, weighted by the transmission factor 
of the path between that element and the receiver. Both the flux and the transmission 
factor are themselves functions of gas temperature. Any one measurement of the nu- 
merical value of this complicated integral cannot provide the local temperature distri- 
bution. In the spectral scanning technique, two or more measurements are made, each 
measure (i.e. , each integral) itself weighted by a function that favors radiation from a 
narrow spectral interval, the interval being different for each measurement. The spec- 
tral intervals are chosen to give greater weight to the radiation from particular portions 
of the optical path. Then, if the additional prior assumption is made that the tempera- 
ture distribution can be characterized by j independent parameters, this distribution 
can be derived, in principle, from j separate and independent measurements. One 
must also know the identity of the gas, the local gas pressure, and the dependence of 
absorption coefficient on temperature, wavelength, and pressure. From dimensional 
considerations, one of the j independent parameters must be a temperature. The other 
parameters may be nondimensional quantities, such as temperature ratios or distance 
ratios, that define the shape of the temperature -distribution curve. 

Each of the "independent measurements" mentioned actually consists of the usual 
set of three flux measurements used in absorption -emission pyrometry: flux from the 
gas alone, flux from a reference source alone, and flux from the reference source after 
it has passed through the gas. 

' The accuracy with which one can determine the temperature distribution will depend 
on the judicious choice of the weighting functions for the individual integrals. Thus, the 
choice of the spectral interval that each measurement favors must maximize the accur- 
acy of temperature measurement. This report addresses itself to the problem of 
spectral -interval selection, for the case j = 2, for profiles that might be encountered in 
turbojet, ramjet, or rocket engines using hydrogen as the principal fuel and oxygen as 
the principal oxidant. 

The problem definition also requires an exact definition of what is meant by "accur- 
acy. " The measure of accuracy that will be used in this report is the error in each of 
the j parameters of the temperature profile. 

The temperature distribution will be assumed to be represented by a curve charac- 
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terized by two parameters - a temperature and a nondimensional variable that defines 
the shape of the curve. The fractional error in temperature due to a given fractional 
error in flux measurement will then be computed for various combinations of wavelength 
and absorption coefficient, in order to identify the combination that yields the smallest 
fractional error in temperature. In order to provide specific numerical results, the 
example of infrared radiation from water, as developed in hydrogen-oxygen or hydrogen- 
air combustion, will be treated in some detail. 


SYMBOLS 

A spectrometer constant, eq. (4) 

b average line spacing 

Cg Planck's second radiation constant 

E error ratio 

F,f factors in error formulas, eqs. (22) and (23) 

G integral of slit function 

g slit function 

H temperature coefficient of absorption coefficient 

k absorption coefficient 

kg absorption coefficient at reference pressure Pq 

L path length 

m number of zones 

spectral radiance 

N. zone-averaged spectral radiance 

n profile shape parameter 

P partial pressure of emitting gas 

P partial pressure of gas z 

z 

T temperature 

T zone -averaged temperature 

T m temperature parameter 

u optical depth 

V,v factors in error formulas, eqs. (20), (24), and (25) 
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X 


distance along path 
y function of x/L 

/3 fine -structure coefficient 

y line half -width 

e absorption coefficient error ratio 

T) " exponent in formula for optical depth 

X ^ wavelength 
Ax effective slit width 

| variable of integration 

to a ^ summations 
r transmission factor 

(p flux 

Subscripts : 

a,b set of measurements 

av average 

L value at x = L 

0 reference 

Superscripts: 
c calculated 

m measured 

ANALYSIS 

Method of Treatment 

In absorption-emission pyrometry , three flux measurements <p^, (p 2 , and are 
made. These represent, respectively, the fluxes from the gas column alone, from the 
source alone, and from the source after passing through the gas column. It is conven- 
ient for computational purposes also to introduce an intermediate variable, the overall 
transmission factor r L = 

In the spectral -scanning method, several sets of such measurements are made at 
different wavelengths and perhaps at different effective slit widths. The number of sets 
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must be not less than the number of independent parameters that characterize the tem- 
perature distribution along the gas column. When the number of sets exceeds the num- 
ber of parameters, the redundancy helps to reduce the statistical inaccuracy of the re- 
sults. In this report, the temperature distribution over the path 0 ^ x :S L will be 
assumed to be characterized by two parameters: a temperature T m and a nondimen - 
sional variable n that characterizes the shape of the temperature profile. Thus, 

T = T(x,T m ,n) (1) 

It will be assumed that the two sets of flux measurements are made, assigned sub- 
scripts a and b, respectively, and that for each set one may select a mean wavelength 
X and an effective slit width AX. The formal error analysis problem is to determine 
the error in the deduced values of T m and n due to errors in the six flux measure- 
ments. The inverse experiment -design problem. is to determine values of X , X^, AX a , 
and AX^ that will lead to low values of the errors | 6T /T m | and | 6n/n | when all 
six flux measurements have the same fractional error \6cp/cp | . The solution of the 
problem requires assumption of a physical model of the mechanisms of spectral emis- 
sion and absorption at each of the two spectrometer settings. Specifically, in the present 
treatment, it requires assumptions of the magnitude and temperature coefficient of the 
absorption coefficient. 

A convenient approach to the problem is to compute the fluxes for an assumed tem- 
perature distribution so that the <p a 's and (p^'s are expressed as functions of T m 
and n. Differentiation then relates the error in each w to errors in T_ and n 
through formulas of the general form 


d cp 




+ 



( 2 ) 


Solution of the simultaneous relations then yields the inverse expression of dT m and dn 
in terms of the errors in the individual cp's. The algebraic expressions are then con- 
verted to a statistical summation of errors by summing the squares of the individual, in- 
dependent error terms. There will be six such terms, corresponding to the six flux 
measurements. Finally, the individual values of fractional error in flux measurements 
are set equal to each other to establish the error ratios 
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( 3 ) 


E(T ) = 
' m 


6T /T 

m 

5 (p/(p 






E(n) = 


6n/n 

6 <p/(p 


General Equation for Gas Radiation 

The flux from a nonisothermal gas column of length L, measured by a receiver 
located at x = 0, is 


*1 


= A 


r r L 

rx 

1 sM 

I N x (X,x)k(X,x) exp - 

/ k(M)d? 

•'A 

Jo t 

Jo *- 

0 


dx dx 


(4) 


where A is a constant of the spectrometer, equal to the product of spectrometer trans- 
mission factor, entrance-slit area, and dispersing-element area (or collimating-element 
area if there is one between the slit and the dispersor), divided by the square of the dis- 
tance between slit and dispersor or collimating element. 

The spectral radiance is given adequately for error analysis by Wien’s approx- 
imation to Planck's equation. The exponential is the gas transmission factor t(X,x). 
Plass (ref. 2) shows that the product k(X,x)-r(X,x) may be written as -3t(X,x)/9x. 

The spectrometer slit function g(X) is appreciably different from zero over only a 
narrow span of wavelength. We define the effective slit width AX as representing this 
span, and represent the integral by G, so that 

G = f g(X)dX = f g(X)dX (5) 

Z J o •'ax 

Equation (4) may be integrated over the wavelength by use of two approximations: 

(1) when AX is small, may be considered independent of X over the span AX, and 

(2) when a spectral band model is available that allows one to effectively integrate 
t(X,x) over the span AX, t(X,x) may be replaced by 
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( 6 ) 


f g(X)r(X,x)dX 
t(x) = ^ 



If the band model is to be valid, AX must include enough spectral lines to produce an 
average r that is not sensitive to a further increase in AX. On the other hand, AX 
must be small enough so that it lies within the band and so that it does not exceed the 
limits within which conditions (1) and (2) are reasonably valid. 

With these approximations, equation (4) becomes 



In the absorption-emission method, where measurements are made of cp-, and of 
t(L) = <p^/<P 2 = t lj a ll three fluxes must be measured with the same value of G. 

Equation (7) cannot be integrated analytically with most temperature profiles found 
in practice, but adequate numerical integration can be performed with any smooth pro- 
file. The path is divided into m isothermal zones, so that equation (7) becomes 

m 

*1 - ag Z ”s (T i ■ T i+i ) (8) 

i=l 

where is the arithmetic average of i and . + j and (i\ - t^) is computed, 
using a band model, as described in the next section and in an example given later on. 


Gas Transmission Function 

The band-model formulation for a nonisothermal gas is based on equation (5) of ref- 
erence 3. This formulation has been justified and applied in reference 4, and summar- 
ized in reference 5. Other band models can be found in the literature (ref. 6). The 
transmission factor r is a function of optical depth u given by 
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where k is the absorption coefficient at temperature T and partial pressure P, and 
j 3 is the fine -structure coefficient. Both k and (3 are functions of T and hence of 
location x along the optical path 0 ^ x < L. The value of 0 may be computed as 

0 = ^221 ( 12 ) 

b 

where y is the line half -width and b is average line spacing. The quantities k, y, 
and b must be known. The quantity 0 av is a mass-weighted average of j3 over the 
entire path. An empirically derived specification of the exponent rj in terms of the 
overall optical depth u L is (ref. 5) 
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Although these equations are used directly for profile determination, approximations 
are necessary to permit presentation of the results of an error analysis in sufficiently 
simple form. The number of independent variables must be limited to those with the 
greatest effect on accuracy. The variables we choose are the overall transmission fac- 
tor t(L) = Tj^ (in preference to the absorption coefficient), the temperature coefficient H 
of the absorption coefficient where H s (T/k) ‘ (3k/9T) , and the wavelength X . The 
wavelength is important because of its strong influence on the spectral radiance, through 
Planck's law. The transmission factor is chosen in preference to the absorption coeffi- 
cient because the former includes the effect of path length. A convenient and adequate 
approximation of H is computed by using the linearization 


k 0< T > = k 0< T m> 




(14) 


where kp is the absorption coefficient at standard reference pressure P Q , and taking 
T = T m - 2(T m - T ), where T av is a temperature that is representative of the major 
flux -contributing portion of the path and T m is a temperature arbitrarily chosen to be 
somewhat higher than T &v . When equation (14) leads to a negative value of k Q (T), 
k Q (T) is to be taken as zero. A reasonable criterion for choosing T av is 

T(x)dx (15) 


To calculate r(x), k g( T m ) is * irst determined, and then equations (14), (10), 
and (9) are used, in that order. 

The value of k o^ T m^ may be °^ taine d from a direct tabulation of this quantity, but 
this method is desirable only when t l is close to unity (a more exact statement of this 
criterion will be given later on). Ordinarily, greater accuracy is achieved by computing 
ko(T m ) from the measurement of t^. The procedure is as follows: 

Given T(x) from an assumed temperature profile, P(T) from the gas composition 
and partial pressure, and from a direct tabulation or from values of y and b, find 
(3 av by equation (11) . From the measured t l s t(L), solve equation (9) for u L and 
then estimate tj from equation (13). Assuming a value of H, equation (14) is used to 
express kp(T) in terms of k g^ T m^ in e( l uation ( 10 )- Integrating equation (10) over the 
entire path length L then yields 
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k 0< T m> = 


/© 

•/n 


U L P 0 


-) 1 + H ( - — -™\ ( dx 


Error Analysis 


Error analysis is effected by differentiation of equation (7) and then replacing dif- 
ferentials of significant quantities by the errors in those quantities. Differentiation of 
equation (7) requires determination of the differentials dN^ and d(dr/dx). The expres- 
sion for dN^, obtained by differentiating Wien's equation, is 

dN x - fe^dT (17) 


where c ^ is Planck’s second radiation constant. 

The quantity dr/dx = r f may be rq>resented as a function of T, k^, and /3 if tabu- 
lated values of k Q are used to determine t(x). Then, one can write 


dT, = 


(18a) 


On the other hand, if the measured is used to establish r(x), one can write 


d T ' = ( — )dT + 


— — — |dTy + ^ d/3 


We shall treat this latter form in detail, and merely summarize, later, the result of 
using equation (18a). 

The last term on the right side of equations (18a) and (18b) will be neglected, for the 
purposes of error analysis, because it is of smaller magnitude than the other two terms. 
We note in passing that 
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but continue to work with for conciseness. 

When equation (7) and its differential are combined with equations (17) and (18) and 
when dT is represented in terms of T m and n through the differential of equation (1) , 
the result may be written as 


dn - f 


dT d <pi 

— = V = v, - — - - v«j dr T 

■m n 2 


( 20 ) 


where f, Vp and V 2 depend on the temperature profile. 

When two sets of measurements are made, identified by respective subscripts a 
and b, there result two equations of the form of equation (20) representing the differen- 
tials in the six flux measurements (through eq. (19)) due to differentials dT m and dn. 
Simultaneous solution of the two equations yields 


dn = 


f V, - f. V 0 
a b b a 

f - f- 
a b 


(21a) 


dT 


m 


m 


V b- V a 

f a- f b 


(21b) 


If all six flux measurements are nowassigned equal fractional error, the corresponding 


errors in T m and n are given by 


where 


E 2 (T ) = — k 
' mr 


F , 2 + F 2 


M 2 


E 2 (n) = 


f 2 F 2 + f 2 F 2 
a b b a 


( f a-' b )= 


F 2 = v 2 t 2 v 2 t 2 


( 22 ) 


(23) 


(24) 
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with subscript a or b assigned to F, r L , f, v^, and Vg. The error ratios are de- 
fined by equation (3). These represent the sensitivity of T m and n to errors in flux 
measurement. The sensitivity is greatest when f ~ f^; hence this difference should be 
held large, and the F-values should be as small as possible. 

Alternative determination of absorption . - In those situat ions when r L is near 
unity, so that the gas is quite transparent, F becomes undesirably large. In that case 
a smaller value of F may be obtained by estimating by a computation based on in- 
trinsic gas properties rather than by obtaining it from experimental flux measurements 
q >2 and cp^- This alternative approach is possible if a tabulation of kg is available and 
if a good a priori estimate of the temperature profile can be made. (This estimate might 
be the profile first determined by the three experimental flux measurements.) Let 
e = (dkQAg)/(d<p/<P) be the ratio of the error of the tabulated kg to the error of flux 
measurement. Then equation (24) can be replaced by 

r2 = v ? + v 2 T L( eInT L) 2 < 25 > 

Accuracy is thereby improved whenever (e In t^J 2 < 2. 


EXAMPLE 

Assumptions and Data for Computation of f and F 

To clarify what magnitudes of error ratios might be achieved in a practical applica- 
tion and what control might be exercised on these error ratios, numerical calculations 
have been made for a specific family of temperature profiles and for the infrared radia- 
tion from HgO. 

The profiles will be assumed to be given by 



where 


y = 




(26) 


(27) 
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and 


T l = 300K, 1000 < T m < 3000K, 0 < n < 1/2 

These profiles are shown in figure 1 . They resemble the transverse velocity and tem- 
perature distributions often found in ducts with axisymmetric flow. 


Profile shape 
parameter, 

V 



Figure 1. - Assumed temperature profiles characterized by two parameters. 


For this profile, dT in equation (17) becomes 


dT = (1 


- y 1/n > dT m + 



- T L )(y l/n 



(28) 


and the values of f, Vj, and Vg in equation (20) become 
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( 29 ) 


<r, = y N. At. 
4 £j i i 



(T. + T. ,) 
T = 1 1+1 

1 2 


AT i = T i - T i + 1 
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where all summations are for i = l to i = m. To compute cr^, and n are held 
constant; to compute a,-, T and n are held constant; to compute u n , r T and T 
are held constant. 

The water radiation will be assumed to lie within the 1.9-, 2.7-, and 6.3- 
micrometer bands. Other assumed limits will be overall transmission factor, 

0.1 < r L < 1; and temperature coefficient of k, -3 £ H < 2. For the profiles of equa- 
tion (26), equation (15) leads to values of T ay that lie between 70 and 100 percent 
of T m . Therefore, a safe rule of thumb for T av is T av «0.85T m . Table I gives 
values of kg and H for wavelengths that are of particular interest, because their use 
can lead to relatively small values of the error ratios. The values in table I are com- 
puted from the tabulations by Ludwig (ref. 7) of values of kgT m /Tg at 1 bar pressure. 
The value of e for the data is on the order of 10. 


TABLE I. - ABSORPTION COEFFICIENT k 0 T m /T 0 AND ITS 
TEMPERATURE COEFFICIENT H FOR WATER AT 
100 KILONEWTONS PER SQUARE METER (1 BAR) 


Wave- 


Temperature parameter, 

T m> K 


length, 

X, 

1000 

2000 

3000 

fim 

k 0 T m /T 0 

H 

k 0 T n/ T 0 

H 

k 0 T m/ T 0 

H 


(a) 

(b) 

(a) 

(b) 

(a) 

(b) 

6.061 

0.237 

-2.04 

0.171 

-0.47 

0.170 

-0.02 

5.970 

.294 

-1.80 

.210 

-.53 

. 202 

-.12 

5.882 

.434 

-.69 

.268 

-1.23 

.173 

-1.51 

5.797 

. 528 

-.62 

.300 

-1.48 

.191 

-1.71 

2.703 

.673 

-.13 

.268 

-2.45 

.150 

-2.36 

2.685 

. 566 

-1.17 

. 194 

-2.60 

.113 

-2.15 

2.597 

. 277 

-2.77 

.161 

-.57 

.154 

-.14 

2.500 

. 186 

.69 

.159 

-. 50 

.161 

.04 

2.484 

. 113 

.75 

.124 

.0 

. 132 

. 18 

2.424 

.012 

1.41 

.035 

1.22 

.060 

1.22 

2.094 

.001 

2.02 

.019 

2.78 

.040 

1.61 

2.020 

.010 

1.99 

.049 

1.42 

.048 

-.09 

1.961 

.053 

1.56 

.048 

-.64 

.040 

-.60 

1.810 

.045 

.56 

.045 

-.11 

.044 

-.10 


a This quantity is tabulated rather than kg to facilitate linear in- 
terpolation between columns. Tp = 273 K; kp = kp(T m ). 

b Evaluated at T = 0.85T . 

dv m 
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The values of y and b are also computed from reference 7. Using values there 
tabulated, one obtains 



The value of b (designated d in ref. 7) can be approximated for HgO by making it inde- 
pendent of wavelength and with a temperature dependence given by 



where 


b Q = 0 . 5 at T q = 1500 K 


For a given amount of fuel and oxidant, the more important partial pressures P 


and can be determined (ref. 8) to allow computation of y 


and, thereby, of 0. 


H„0 

The 


approximations make /3(T) the same at all wavelengths in the range of X . 

The error -ratio parameters f (eq. 29)) and F (eqs. (24) and (29)) are listed in 
tables II and III for a particular wavelength in two of the water bands; interpolation may 
be used for any other wavelength. The quantity f is affected to a lesser degree by 
than it is by X or H. On the other hand, parameter F is more strongly affected by 
t l than it is by X or H. 

To find values of f and F by use of tables II and III , r L must be estimated from 
the optical depth u L - Figure 2 shows the relation between t l and the optical depth u L 
of an isothermal gas column, in accordance with equation (9). The optical depth u^ of 
a nonisothermal gas column is an integral of the local value of k over the path L. For 
the profile of equation (26) , this integral is approximated by 


U L = k 0 (T m ) 




(32) 


for H ^ T m /(T m " TjO. Because a larger value of H would give a negative value of k 
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Temper- 

Profile 

Temper- 

ature 

shape 

ature 

pa ram- 

pa ram- 

coeffi- 

eter , 

eter, 

cient of 

T , 
m’ 

n 

absorp- 

K 


ti on 
coeffi- 
cient, 
H 


Transmission factor at x = L, 


0 . 

1 

0 . 

4 

0 . 

8 


0.95 



Wavelength, X, p 

m 




1.9 

6.3 

1.9 

6.3 

1.9 

Li 

.3 

1.9 

6.3 



Reciprocal wavelength, 

iA 



0. 526 

0. 159 

0.526 

0.159 

0. 526 

0 

.159 

0.526 

0.159 

f 


1000 

0. 1 

-2 

1 . 

5 

0. 

7 

1 . 

9 

1 . 

0 

1 . 

. 5 

0 . 

.6 

1 . 

.6 
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1 . 
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9 

2. 

.4 

1. 

3 

2. 

2 

1 . 

1 

2. 

4 

1.2 



1 

2. 

1 

1 . 

2 

2 

.8 

1 . 

.7 

3. 

.0 

1 . 

6 

3. 

3 

1.8 



2 

2. 

5 

1. 

6 

3 

.5 

2. 

.5 

4. 

.5 

3. 

1 

5. 

0 
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.6 
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.2 

4. 

9 
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.5 

5. 

3 
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TABLE HI. - VALUES OF (1 - t l )F FOR SELECTED VALUES 
OF OTHER PARAMETERS 


[Linear interpolation may be used to obtain intermediate values. 


Temper- 
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T . 
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Profile 

shape 

param- 

eter, 
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1000 

0.1 


0.25 


0.5 

2000 

0.1 


0.25 


0.5 

3000 

0. 1 


0.25 


0.5 


ature 
coeffi- 
cient of 
absorp- 
tion 
coeffi- 
cient, 

H 


Transmission factor at x = L, 
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1.9 6.3 1.9 6.3 


6.3 1.9 6.3 
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Optical depth, u L 

Figure 2. - Overall transmission factor of isothermal gas column. 


at T t , the quantity H should arbitrarily be assigned the limiting value T /(T - T t ) 

■Lj m m l 

if H is actually larger than this value. Pressure P of the emitting gas is a repre- 
sentative value over the temperature range from to T . Representative values 
of P and /3 are their values at the temperature T used in the determination of H. 

3. V 

Procedure for Computing Error Ratios 

The steps in the computation of error ratios are 

(1) Forgiven T m , n, H, X, L, and partial pressure P, obtain values of f from 
table I and then use equation (32) to compute the optical depth u^. 

(2) Compute temperature T = 0. 85 T . 

av m 

(3) Calculate |3 from equation (12) at temperature T ay . 

(4) Obtain from figure 2. 

(5) Obtain f and F from tables II and III, choosing the procedure (measurement or 
calculation) that yields the smaller value of F. For the procedure that uses calculation, 
use the value of F at which = exp(- \Z2/e). 

(6) Calculate the error ratios by equations (22) and (23). 


Computed Error Ratios 


The error ratios E(T ) and E(n) are plotted in figure 3 against the path length L 
for three cases: both measurements in the same band (2.7 pm); measurements in 
adjacent bands (2.7, 6.3 pm and 1.9, 2.7 pm); measurements in widely different bands 


(1.9, 6.3 pm). Values of 
m 


m. 


Tj^ is the ratio of the measured fluxes 


and along the curves are listed in table IV, where 

<p 2 and <p 3 


and C T^ is calculated from a 
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Figure 3. - Error ratios for several wavelength pairs as function of geometric path length. The choice of measured or computed values of was made so as to 
minimize the error ratio. 
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TABLE IV. - TRANSMISSION FACTORS FOR CURVES 


OF FIGURE 3 WITH n = 0. 5 


Tempera- 
ature pa- 
rameter, 
T 

m’ 

K 

Wave- 

length, 

pm 


Path length 

L 


4 

10 

20 

50 

100 

Transmission factor at x 

L 

1000 

5.970 

0.85 

0.74 

0.63 

0.47 

0.33 


2.685 

.79 

.66 

.54 

.37 

. 24 


2.500 

.93 

00 

-q 

.79 

.67 

.55 


2.424 

a . 995 

a . 988 

a . 976 

a . 95 

3 . 91 


1.961 

3 . 979 

a . 95 

3 . 91 

.83 

.75 

2000 

5.797 

a 0. 87 

0.73 

0.58 

0.37 

0.22 


2.685 

3 . 88 

.76 

.62 

.41 

.25 


2.484 

3 . 961 

a .91 

. 84 

.68 

.53 


2.424 

3 . 993 

3 . 982 

3 . 964 

3 . 917 

.85 


2.020 

3 . 990 

3 . 975 

3 . 95 

3 . 89 

.80 

3000 

5.882 

a 0. 93 

0.85 

0.73 

0.50 

0.32 


2.685 

a . 95 

3 . 88 

.78 

.57 

.39 


2.500 

3 . 963 

a . 91 

.84 

.67 

.49 


2.424 

3 . 991 

a .980 

3 . 960 

3 . 91 

.83 


2.094 

a . 995 

a .986 

3 . 973 

3 . 94 

00 

CO 


a This quantity is c r L ; others are 


tabulation of absorption coefficients. The error ratio E(T m ) is smaller at lower temper- 
atures because (1) |f a - f ^ | is greater, and (2) F decreases with decreasing r L , which 
decreases with decreasing temperature, unless path length also decreases. For ex- 
ample, at a path length of 10 centimeters the lowest obtainable value of E(T m ) is about 
6 at 3000 K and about 2 at 1000 K. The lowest obtainable value of E(n) is relatively inde- 
pendent of temperature. 


PROCEDURE FOR SELECTING OPTIMUM WAVELENGTHS 


Given a gas of path length L, pressures 


p o 5P and 



wall temperature 


T 


L’ 


and prior estimates of T m and n, the following steps may be used to select the pair of 
wavelengths X a> x^ that will minimize the error ratios E(T m ) and E(n): 

(1) Use table II to select approximate values of X a , X^, and corresponding param- 
eters H to minimize the error ratio in accordance with the selection rule to keep 
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|f - f^| large, and in accordance with considerations of convenience such as availa- 
bility of detectors, filters, or monochromators. Use step (2) to confirm that the chosen 
combination of H and X exists. 

(2) By use of table I, choose wavelengths X , X^, at which the desired values of H 
occur. The error ratio is minimized when the longer wavelength has the more negative 
value of H, and the larger value of k. 

(3) Compute error ratios by the steps listed under ’’Procedure for Computing Error 
Ratios . ' ' 


EFFECT OF LARGE ERROR RATIO ON PROFILE COMPUTATION 

In application of the spectral -scanning method, a profile is computed using flux mea- 
surements and an iterative procedure. The Newton-Raphson iteration method works 
well. When measured flux cp ^ and transmission factor were assumed to be in 
error, trial examples showed convergence of the iteration in a small number of steps 
when the error ratios were small, but the iteration did not always converge when the 
error ratios were large. Instead, a stable oscillation existed between two sets of values 
of T m and n as an approximate solution. The oscillation could be eliminated (ref. 9) 
by use of the damped least squares method (ref. 10). However, this method was not 
necessary when the error ratio was small, as required for practical use of the spectral 
scanning method. 

The stability of the iterative computational procedure depends on (a) the complexity 
of the assumed profile shape, that is, on the number of independent parameters assumed 
to characterize the profile and (b) the magnitude of the error ratios (error in parameter 
divided by error in flux measurement). For a given measurement error S<p/cp, the sta- 
bility of the iterative solution of the simultaneous equations becomes less likely as the 
assumed profile becomes more complex. For a given assumed profile shape, the solu- 
tion of the simultaneous equations resulting from the presumably independent flux fnea- 
surements is more likely to be unstable when the measurement errors are larger. 

The present error analysis, with an example of a constant -pressure gas stream, has 
shown that two measurements in the infrared must be made at carefully selected wave- 
lengths to avoid excessive random error in the determination of the profile. Other cal- 
culations, which will not be detailed here, show that three measurements and an assumed 
three -parameter profile produced a much greater random error. In general, with both 
two- and three -para meter profiles, the random error is reduced and the prospects for a 
stable solution are increased, if measurements are made in both the infrared and the 
ultraviolet. 

When spectral -scanning was applied to measure atmospheric temperature profiles 
from satellites, (refs. 11 and 12), an analysis showed that nominally independent flux 
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measurements may have a high degree of interdependence. In this application, there 
were two favorable circumstances: (1) The probable profile was known to a good approxi- 
mation from accumulated statistical data, and (2) the pressure gradient along the line of 
sight, from satellite to Earth, was more conducive to accuracy than the constant- 
pressure condition assumed in this report. Nevertheless, a maximum of no more than 
six independent measurements was found feasible for that application. A larger number 
caused mathematical instability in the computations, and actually increased the random 
error. 


CONCLUDING REMARKS 

The preceding discussion has treated errors resulting from uncertainty in flux mea- 
surement. Additional sources of error lie in the uncertainty of knowledge of the char- 
acteristics of the gas itself; these include the validity of the band model and, through it, 
the accuracy of the band-model parameters k and j3. These are affected also by un- 
certainties in knowledge of gas composition. On the other hand, the two wavelengths and 
the corresponding Planck function are known with negligible error. The use of widely 
separated wavelengths generally results in a desirably larger value of the quantity 
l f a " %l ■ Such a large value of the difference is less dependent on the value of H or 
t l or on uncertainty in the band model computation of r(x) by equation (9) . Large error 
ratios may also lead to slow convergence of the iterative procedures ordinarily used in 
application of the spectral -scanning method. 

The examples show that E(T m ) increases as T m increases and is quite large at 
short path lengths. The resulting error ratio may be excessive. It is therefore, appro- 
priate to consider more desirable alternatives. One such alternative is to dispense with 
the restriction to infrared radiation and to examine the possible advantages of using 
ultraviolet radiation as well as infrared. At temperatures above 2000 K the pressure of 
OH is sufficient to permit the measurement of its radiation in the ultraviolet. It is then 
possible to make one temperature measurement with OH radiation and only one temper- 
ature measurement with HgO radiation. The wavelength separation between ultraviolet 
and infrared is large. This wavelength separation results in a large difference of param- 
eters |f & - f a | • The short wavelegnth of OH reduces the value of F as indicated by 
table HI . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 17, 1972, 

764-74. 
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